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Abstract 

Here, we peruse cosmological usage of the most promising candidates of dark energy in the 
framework of f{R) theory. We reconstruct the different f{R) modified gravity models in the 
spatially flat FRW universe according to the ordinary and entropy-corrected versions of the 
holographic and new agegraphic dark energy models, which describe accelerated expansion 
of the universe. We also obtain the equation of state parameter of the corresponding f{R)- 
gravity models. We conclude that the holographic and new agegraphic /(i?)-gravity models 
can behave like phantom or quintessence models. Whereas the equation of state parameter 
of the entropy-corrected models can transit from quintessence state to phantom regime as 
indicated by recent observations. 
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1 Introduction 



Astrophysical data from type la supernovae, cosmic microwave background radiation and large 
scale structure have provided convincing evidence for the present observable universe to be in 
the phase of accelerated expansion [1]. One explanation for the cosmic acceleration is the dark 
energy (DE), an exotic energy with negative pressure. The origin and nature of DE is still a 
mystery. A great variety of DE models have been proposed (for review see [2, 3]), but most of 
them are not able to explain all features of the universe, or are artificially constructed in the 
sense that it introduces too many free parameters to be able to fit with the experimental data. 

The holographic DE (HDE) is one of interesting DE candidates which was proposed based 
on the holographic principle [4, 5, 6]. It was shown in [7] that in quantum field theory, the UV 
cut-off A should be related to the IR cut-off L due to limit set by forming a black hole. Following 
this line, Li [8] argued that the total energy in a region of size L should not exceed the mass of 
a black hole of the same size, thus L^pA < LMp, where p\ is the qiiantum zero point energy 
density caused by UV cut-off A and Mp is the reduced Planck Mass Mp^ = SirG. Also Li [8] 
showed that the cosmic coincidence problem can be resolved by inflation in the HDE model, 
provided the minimal number of e-foldings [8]. The HDE models have been studied widely in 
the literature [9, 10, 11, 12, 13]. Obviously, in the derivation of HDE, the black hole entropy 
Sbh plays an important role. As is well known, usually, Sbr = A/{4G), where ^4 ~ is the 
area of horizon. However, in the literature, this entropy-area relation can be modified to [14] 

^BH = ^ + alnA + ^, (1) 

where a and /3 are dimensionless constants of order unity. These corrections can appear in the 
black hole entropy in loop quantum gravity (LQG) [15]. They can also be due to thermal equi- 
librium fluctuation, quantum fluctuation, or mass and charge fluctuations (for a good review see 
[15] and references therein). More recently, motivated by the corrected entropy-area relation (1) 
in the setup of LQG, the energy density of the entropy-corrected HDE (ECHDE) was proposed 
by Wei [15]. 

Recently, the original agegraphic dark energy (ADE) model was proposed by Cai [16]. The 
ADE model assumes that the observed DE comes from the spacetime and matter field fluctua- 
tions in the universe [16]. Following the line of quantum fluctuations of spacetime, Karolyhazy 
et al. [17] discussed that the distance t in Minkowski spacetime cannot be known to a better 

accuracy than St ~ tp'^i^/"^ where tp is the reduced Planck time. Based on Karolyhazy relation, 
Maziashvili [18] discussed that the energy density of metric fluctuations of the Minkowski space- 
time is given by p\ ~ l/(tpt^) Mp/t'^. Based on Karolyhazy relation [17] and Maziashvili 
arguments [18], Cai proposed the original ADE model to explain the accelerated expansion 
of the universe [16]. Since the original ADE model suffers from the difficulty to describe the 
matter-dominated epoch, the new ADE (NADE) model was proposed by Wei & Cai [19], while 
the time scale was chosen to be the conformal time instead of the age of the universe. It was 
found that the coincidence problem could be solved naturally in the NADE model [20]. The 
ADE models have given rise to a lot of interest recently and have been examined and studied 
in ample detail in [21, 22, 23, 24]. More recently, very similar to the ECHDE model, the energy 
density of the entropy-corrected NADE (ECNADE) was proposed by Wei [15] and investigated 
in ample detail by [25]. 

Another explanation for the cosmic acceleration and what we get from observational data is 
the different approach, so called "modified gravity" . This scheme possesses the relevant feature 
that experimental data turn out to be naturally interpreted without the need of additional 
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components like DE (for review see [26]). In the situation when Einstein general relativity 
(GR) cannot naturally describe the DE epoch of the universe the search of alternative, modified 
gravity which is consistent with solar system tests/observational data is of primary interest 
[27]. A very popular modified gravity model is the so-called /(i?)-gravity [28, 29, 30] where 
f{R) is an arbitrary function of the scalar curvature R. Actually, among extended theories of 
modified gravity, /(i?)-gravity represents a viable alternative to DE and naturally gives rise to 
accelerating singularity-free solution in early and late cosmic epochs [31]. For instance, the form 
f{R) = R + R? + 1/R can predict the unification of the early-time inflation and late-time cosmic 
acceleration in the standard metric formulation [32]. 

Viewing the modified /(i?)-gravity model as an effective description of the underlying theory 
of DE, and considering the ordinary and entropy-corrected versions of the HDE and NADE 
scenarios as pointing in the direction of the underlying theory of DE, it is interesting to study 
how the /(i?)-gravity can describe the HDE, ECHDE, NADE and ECNADE densities as effective 
theories of DE models. This motivated us to establish the different models of /(i?)-gravity 
according to the ordinary and entropy-corrected versions of the HDE and NADE scenarios. This 
paper is organized as follows. In section 2, we review the theory of /(i?)-gravity in the metric 
formalism. In sections 3 to 7 we reconstruct the different /(-R)-gravity models corresponding 
to the HDE, ECHDE, NADE and ECNADE models, respectively. Section 8 is devoted to our 
conclusions. 



2 The theory of f{R) modified gravity 



The general /(/2)-gravity action is given by [27, 28] 

-R + f{R) 



S 



(2) 



2]^2 ' ^matter 

where = Mp"^ = SnG and h = c=l. Also G, g, R and 

-^matter SJ"© the gravitational constant, 
the determinant of metric g^,^, the Ricci scalar and the lagrangian density of the matter inside 
the universe, respectively. 

One notes that there are in fact two strategies in f{R) theories: the metric formalism, where 
the action is varied with respect to the metric only; and the Palatini formalism [33], where the 
metric and the connection are treated as two independent variables with respect to which the 
action is varied [34]. It is only in Einstein gravity f{R) = that both approaches reach the 
same result. In general f{R) theories, the problem of which approach should be used is still an 
open question and the final solution may be determined by further observations and theoretical 
development [34]. Here like [34] we use the metric formalism. 

Taking the variation of the action (2) with respect to the metric g^i^, one can obtain the field 
equations as [27, 28] 

R,,-^Rg,. = e{T(^^+T(^^), (3) 

where 

k^T(^^ = Ig^ufiR) - R^.f'{R) + (V^V. - g^.n)f{R). (4) 

Here R^y and T^™^ are the Ricci tensor and the energy-momentum tensor of the matter, respec- 
tively. Also the prime denotes a derivative with respect to R. 

Now if we consider the spatially fiat FRW metric for the universe as 

3 

= _dt2 + a2(i)^(da;^)2, (5) 
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and taking Tv^^^ = diag{—pm,Pm,PrmPm) for the energy-momentum tensor of the matter in the 
prefect fluid form, then the set of field equations (3) reduce to the modified Priedmann equations 
in the framework of /(i?)-gravity as 



where 



= Pm + PR, 



fe2 



(2H + 3if^ 



-{Pm+PR), 



--/(i?) + 3{H + H^)f\R) - 18{iH^H + HH)f"{R) 



(6) 
(7) 

(8) 



PR 



A;2 



and 



■-fiR)-{H + 3H')f{R) 

+6{8H^H + 6HH + AH^ + H)f"{R) + S6{H + AHHff\R) 



R = 6{H + 2H^). 



(9) 



(10) 



Here H = a/a is the Hubble parameter and the dot denotes a derivative with respect to cosmic 
time t. Also pR and pR are the curvature contribution to the energy density and pressure. 
The energy conservation laws are still given by 



Prn + SH{prn + Pm) = 0, 
PR + 3H{pR+PR)=Q. 



[ID 



(12) 



In the case of f{R) = 0, from Eqs. (8) and (9) wc have /9r = and pR = 0. Therefore Eqs. (6) 
and (7) transform to the usual Friedmann equations in GR. 

The equation of state (EoS) parameter due to the curvature contribution is defined as [35] 



PR 

OJR= — 

PR 



Hf'iR) + 3{3HH - AH^H + iH^ + H)f"{R) + 18{H + AHHff'iR) 



f{R) - 6{H + H^)f'iR) + 36{4:H^H + HH)f"{R) 

Note that for a f{R) dominated universe, Eq. (6) yields 

3 o 



(13) 



(14) 



Taking time derivative of the above relation and using the continuity equation (12), one can get 
the EoS parameter as 

2H 

which shows that for the phantom, ur < — 1, and quintessence, ujr > — 1, dominated universe, 
we need to have H > and H < 0, respectively. 

For a given a = a{t), by the help of Eqs. (8) and (9) one can reconstruct the /(i?)-gravity 
according to any DE model given by the EoS pR = Pr{pr) or pR = pR{a). There are two classes 
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of scale factors which usually people consider them for describing the accelerating universe in 
f{R), f{g) and f{R,g) modified gravities [36]. 
The first class of scale factor is given by [36, 37] 

a{t) = ao{ts - t<ts, h>0. (16) 

Using Eqs. (10) and (16) one can obtain 

h 



H " 



ts-t [6(2/1 + 1 



-R 



1/2 

, H = H^/h, (17) 



which H = H'^ /h > shows that the model (16) is correspondence to a phantom dominated 
universe. This is why that in the literature the model (16) is usually so-called the phantom scale 
factor. 

For the second class of scale factor defined as [36] 

a{t) = aot^, h>0, (18) 

one can obtain 

h 



t 



-R 



1/2 

, H = -H^/h, (19) 



_6{2h - 1) 

which H = —H'^/h < reveals that the model (18) describes a quintessence dominated universe. 
Hence, the model (18) is so-called the quintessence scale factor in the literature. 

In sections 3 to 6 using the two classes of scale factors (16) and (18), we reconstruct the 
different /(i?)-gravities according to the HDE, ECHDE, NADE and ECNADE models. 

3 f{R) reconstruction from HDE model 

Here we reconstruct the /(i?)-gravity according to the HDE scenario. Following Li [8] the HDE 
density in a spatially flat universe is given by 

where c is a numerical constant. Recent observational data, which have been used to constrain 
the HDE model, show that for the flat universe c = 0.818^0 097 [38]. Also Rh is the future event 
horizon defined as 

/•°° dt da 

For the first class of scale factor (16) and using Eq. (17), the future event horizon R^ yields 



dt ts-t 1 Qhi2h + l) 
Replacing Eq. (22) into (20) one can get 

Substituting Eq. (23) in the differential equation (8), i.e. pR = pA, gives the following solution 

f{R) = X+R"'+ + X-R""- + 7c R, (24) 
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where 

3 + h±Vh'^ -lOh + 1 
m± = ^ , (25) 

and 

7.^-^. (26) 

Also X± are the integration constants that can be determined from the necessary boundary 
conditions. Fohowing [39] the accelerating expansion in the present universe could be generated, 
if one consider that f{R) could be a small constant at present universe, that is 

f{Ro) = -2Ro, (27) 

/'(i?o) ~ 0, (28) 

where Rq ~ (10~^^eV)^ is the current curvature. Applying the above boundary conditions to 
the solution (24) one can obtain 



7c(Tn_ — 1) + 2m- 
(m_|_ — m-)R^'^ 



A+ = (Lt:; . (29) 



A_ = 7e(rn+-l) + 2m+ _ 

Replacing Eq. (24) into (13) and using (17) one can get the EoS parameter of the holographic 
/(i?)-gravity model as 

0JR = -l-^, h>0, (31) 

which corresponds to a phantom accelerating universe, i.e. lor < — 1. Recent observational data 
indicates that the EoS parameter ojr at the present lies in a narrow strip around lor = — 1 and 
is quite consistent with being below this value [3]. 

For the second class of scale factor (18), using Eq. (19) the future event horizon R^ reduces 

to 

f°°dt t 1 /6/i(2/i- 1) ^ ^ 

where the condition /i > 1 is obtained due to have a finite future event horizon. If wc repeat 
the above calculations then one can obtain the both of f{R) and w/j corresponding to the HDE 
for the second class of scale factor (18). The result for f{R) is same as (24) where now 

3-h±VWTWhTl 
m± = , (33) 

7e = -^^. (34) 

Also the EoS parameter is obtained as 

ujR = -l + ^, h>l, (35) 

which describes an accelerating universe with the quintessence EoS parameter, i.e. — 1 < ur < 
-1/3. 
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4 f{R) reconstruction from ECHDE model 



Using the corrected entropy-area relation (1), the energy density of the ECHDE can be obtained 
as [15] 



PA 



n4 ' 



(36) 



where a and ^ are dimensionless constants of order unity. In the special case a = ^ = 0, the 
above equation yields the well-known HDE density (20). Since the last two terms in Eq. (36) 
can be comparable to the first term only when is very small, the corrections make sense only 
at the early stage of the universe. When the universe becomes large, ECHDE reduces to the 
ordinary HDE [15]. 

For the first class of scale factor (16), substituting Eq. (22) into (36) yields 



PA 



:R + 



{h + l) 



2A;2/i(2/i + 1) 36/i2(2/i + l)2 



/ Qh{2h + 1) 



(37) 



Solving the differential equation (8) for the energy density (37) gives the entropy-corrected 
holographic /(i?)-gravity as 



f{R) 



A;2(/i + l)4 



Mh?{2h + l) 



a 



//i-5\_ f Qh{2h + l) 



-I5\R\ 



(38) 



where m± and 7c are given by Eqs. (25) and (26), respectively. Also \± are determined from 
the boundary conditions (27) and (28) as 



A. 



1) + 2m. 



{m+-m.)Rl 



,m+ — 1 



+ 



l)4(m_ - 2) 



54/i2(2/i + l)(m+ -m_)i?^+ ^ 
h-b 



a 



1 / 6h{2h + l) \ 

m_-2 "^[k^ih + iyRoJ 



(39) 



A_ 



--2 



X 



7c(m+ -1) + 2m+ P{h + l)^(m+ - 2) 

(m_ - m+)i?^-~^ Uh^{2h + l)(m_ - m+)i?^ 

' 1 / 6h{2h + l] 



In 



(40) 



3 m+-2 \k'^{h + iyRo. _ 

Substituting Eq. (38) into (13) one can get the EoS parameter of the entropy-corrected holo- 
graphic /(i?)-gravity model as 



= -1-:^ 1 + 



fc2(/i + l)2 




+ /3 


R 


18c2/i(2/l + l) + it2(/i+l)2 


a In 


( 6h{2h+l) \ 


+ P 


R 



h>0, (41) 



which can be rewritten using the first relation of Eq. (17) as 



= -1 - 



-a + 2a In 



1 3c2 (^^) % 2a In (,^) + P 



+ /3 



h>0. 



(42) 
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The above relation shows that the EoS parameter is time-dependent and in contrast with con- 
stant EoS parameter (31), it can justify the transition from quintessence state, cor > —1, to 
the phantom regime, cj/j < —1, as indicated by recent observations [40]. Note that if we set 
a = /3 = then Eqs. (38) and (42) reduce to (24) and (31), respectivefy. 
For the second class of scale factor (18), the result of f{R) is obtained as 



f{R) = \+R"'+ + X.R"'- + 7ci? 
A:2(/i- 1)4 



where 



+ 



54/i2(2/i- 1) 



a 



h + 5\ , /6/i(2/i-l) 
+ m ' 



W{h-l)^R 



A4- = 



1) + 2m_ 



+ 



k'^{h-l)\m--2) 



(m+ - m-)Ro+~^ 54/i2(2/i - l)(m+ - m_)i?^+~^ 



a 



h + 5 



+ 



In 



m- 



/ 6h{2h ~ 1) 
Wih-lfRo 



+ (3\R\ 



+ /5 



(43) 



(44) 



+ 



eih-l)\m+-2) 



7c(m+ - 1) + 2m-| 

(m_ - m+)i?^-~^ ' 54/i2(2/i - l)(m_ - m+)i?^-~^ 



X 



a 



h + 5 



+ 



1 I . fjhi^h^ 
m+-2^ \k^{h-l)^Ro 



(45) 



and the parameters m± and 7c are given by Eqs. (33) and (34), respectively. Also the EoS 
parameter is obtained as 



-1 + — < 1 + 
3h^ 



-« + 2«ln(^(^)+/3 



h>l. 



(46) 



Here also to have a finite Rh, the parameter h should be in the range oi h > 1. One notes that 
the EoS parameter (46) is also dynamical and in contrast with constant EoS parameter (35), it 
can accommodate the transition from lor > — 1 to lor < — 1 at recent stage. 



5 f{R) reconstruction from NADE model 

Following [19], the energy density of the NADE is given by 



PA 



3n^ 



(47) 



where the numerical factor Sri^ is introduced to parameterize some uncertainties, such as the 
species of quantum fields in the universe, the effect of curved spacctimc (since the energy density 
is derived for Minkowski spacetime), and so on. It was found that the coincidence problem could 
be solved naturally in the NADE model provided that the single model parameter n is of order 
unity [20]. The joint analysis of the astronomical data for the NADE model in flat universe 
gives the best-fit value (with la uncertainty) n = 2.716I!Iq;JJ9 [20]. Also 77 is the conformal time 
of FRW universe, and given by 

dt f da 



(48) 
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For the first class of scale factor (16), the conformal time r] by the help of Eq. (17) yields 

dt {ts - 1 



r? 



R 



(49) 



PA = ,:iz: (50) 



It a ao(^ + l) ao(/i + l) 
Substituting Eq. (49) into (47) one can obtain 

A;2(6/i(2/i + 1))' 

Solving the differential equation (8) for the energy density (50) reduces to 

f{R) = X+R"'+ + X.R""- + 7n R''+\ (51) 

where 



^ = + (52) 

h'^{h + 2){6h{2h + l)f' 

and m± are given by Eq. (25). Using the boundary conditions (27) and (28) the parameters X± 
are determined as 

^ 7n(m- -ih + 1))R^ + 2m- ^^^^ 
(m+ — m_)i?^+~^ 
^_ ^ 7n(m+ -{h + 1))R^ + 2m+ 

/ N. j-,m- — l ' ^ ' 

(m_ — m-^.)RQ 

Replacing Eq. (51) into (13) one can get the EoS parameter of the new agegraphic /(i?)-gravity 
model as 

coR = -l-^^^, h>0, (55) 

which like the EoS parameter of the holographic /(i?)-gravity model (31), it always crosses the 
phantom-divide line, i.e. lur < — 1. 

For the second class of scale factor (18), using (19) the conformal time r] is obtained as 



dt t^-^ 1 



r]= - 



Qh{2h - 1) 



R 



\<h<l, (56) 



/o a ao(l - h) ao(l - h) 

where the condition ^ < /i < 1 is necessary due to have a real finite conformal time. The result 
of f{R) is 

f{R) = X+R"'+ + X-R""- + jn R^~^, (57) 

where 

7n(m- + h- l)i?o ^ + 2m_ 
+ } N ' 1^8 j 

^ 7n(m+ + /t - l)i?o ^ + 2m+ 
(m_ — m_|_ j/tQ 



with 

n^ag(/i-l) ^ 
/i2(/i-2)(6/i(2/i- 1)) 
and the parameters mj- are given by Eq. (33). 
Also the EoS parameter is obtained as 



which shows a quintessence-like EoS parameter with — 1 < ojr < —1/3. 
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6 f{R) reconstruction from ECNADE model 



With the help of quantum corrections to the entropy-area relation (1) in the setup of LQG, the 
energy density of the ECNADE is given by [15] 



PA 




(62) 



which closely mimics to that of ECHDE density (36) and is replaced with the conformal 
time r). Here a and P are dimensionless constants of order unity. In the special case a = P = 0, 
the above equation yields the well-known NADE density (47). 

For the first class of scale factor (16), substituting Eq. (49) into (62) one can get 



PA 



+ 



(6/i(2/i + 1)) 



2/1+2 



h+l 



a In 



{6h{2h + l)) 



+ /3 



(63) 



Solving the differential equation (8) for the energy density (63) yields 



f{R) = X+R"'+ + X.R""- + 7„ R''+^ 



3/i(3 + lOh + 6/i2)(6/i(2/i + 1)) 



2/t+l 



a 



(7/i + 5)(/i + l) 



h+l 



+ In 



(6/i(2/i + 1)) 
A;2a2(/i + l)2i?fc+i 



+ /3 R"'+\ 



(64) 



3 + lO/i + 6/i2 

where m± and 7„ are given by Eqs. (25) and (52), respectively. Also X± are determined from 
the boundary conditions (27) and (28) as 

7„(m_ - {h + l))R^ + 2m- 



+ 



(m+ - m.)R'^^-' 

e4{2h + 2 - m_)(/i + l)^Rf^^-"'+ 



3/i(3 + lO/i + 6/i2)(6/i(2/i + 1)) 



2fe+l, 



a 



(7/i + 5)(/i + l) 



+ 



h + l 



A- 



3 + 10/i + 6/i2 ■ m_-2-2/i 

7n(m+ - (/i + l))fi(^ + 2m 
(m_ — m+)/tg 



-Fin 



^ (6fe(2fe + 1))^+^ \ 
^Pa2(h + l)24+ijj 



+ /3 



(65) 



+- 



e4{2h + 2-m+){h + l)4i?f +2-"^- 



3/i(3 Wh + 6/i2)(6/i(2/i + 1)) 



2h+l, 



(7/i + 5)(/i + l) 



+ 



h + l 



+ ln 



(6/i(2/i + 1)) 



/i+i 



+ /3 



(66) 



3 + 10/i + 6/i2 ' m+-2-2/i '"^\k^al{h + lfR'f^+\ 

Replacing Eq. (64) into (13) yields the EoS parameter of the entropy-corrected new agegraphic 
/(i?)-gravity model as 



UR = -1- 



2{h + l) 
3h 



—a -|- a In 



1 + 



(6/i(2fe+l))''^^ 
fc2ag(/i+l)2flh+i 



+ /3 



3n^ 



fe2ag(/i+l)i 



6/t(2/t+l)xh+l 



-|- a In 



(6/t(2/i-l)) 
fe2a^(/i+l)2il'»+i 



/i > 0. (67) 
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It can be rewritten using the first relation of Eq. (17) as 



2(/t + l) 
3h 



1 + 



> , h>0, (68) 



which is time-dependent and in contrast with constant EoS parameter (55), it can justify the 
transition from ujn > —1 to cjr < —1. Note that if we set a = /3 = then Eqs. (64) and (68) 
reduce to (51) and (55), respectively. 

For the second class of scale factor (18), the result of f{R) is obtained as 

ea^{h - 1)4 



f{R) = + + 7„ R^-^ + 



a 



3/i(3 - lO/i + 6/i2)(6/i(2/i - 1)) 

l-h ' 



l-2h 



where 



3 - lO/i + 6/i2 
7n(m- + {h- l))Ro^ + 2m- 



(7/i - - 1) ^ / (6/i(2/i - 1)) 



+ /3}r: 



2-2h 



(69) 



4 jj2— 2/i— m+ 



3/i(3-10/i + 6/i2)(6/i(2/i- 1))^ ^''(m+-m 



a 



{7h-5){h-l) 



+ 



1 - h 



l-h 



+ In 



A_ 



3-10/i + 6/i2 m_-2 + 2/i 
7n(m+ + (/t - l))it:o ^ + 2m+ 

/ N. j-,m- — l 

(m_ — m-\-)RQ 
3h{3 - lOh + Qh?){ioh{2h - \)f~'^^{m- - m+) 



{<dh{2h - 1)) 
k^al{h-lfRl-\ 



+ /3 



(70) 



a 



+ 



1-^ 



(7/t-5)(h- 1) 
3 - lO/i + 6/i2 ' m+-2 + 2/i 



+ ln 



' (6/i(2/i-l)) 
^Pa§(/.-l)2i2i-\ 



(71) 



and the parameters m± and 7„ are given by Eqs. (33) and (60), respectively. Also the EoS 
parameter is obtained as 



-, 2(1 -/i) , ^ 



-« + 2aln(,^(A)i-'^)+^ 



3nK^^(f)^-') +2«ln(,^^(A)i-^)+/?^ 



-</.<l. 

(72) 

Here also to have a real finite conformal time 77, the parameter h should be in the range of 
I < < 1. Contrary to the constant EoS parameter (61), the dynamical EoS parameter (72) 
can accommodate the transition from ur > — 1 to ojr < — 1 at recent stage. 

7 f{R) reconstruction in de Sitter space 

The scale factor in de Sitter space is defined as 

a{t) = aoe^*, H = constant, (73) 
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which can describe the early-time inflation of the universe [36]. Using Eqs. (10) and (73) one 
can obtain 



Then Eq. (8) takes the form 

4k^PR = -2f{R) + Rf'{R). (75) 

Also the EoS parameter (13) yields lor = —1 which behaves like the cosmological constant. 
For the scale factor (73), using Eq. (74) the future event horizon reduces to 



Substituting Eq. (76) into HDE density (20) one can get 

PA = ^R. (77) 

Replacing Eq. (77) in the differential equation (75), i.e. pR = pA, gives the solution 

f{R) = \R^ - c^R, (78) 

where A is an integration constant. Note that in order to generate the inflation at the early 
universe as in Starobinsky's model [36, 41], one may require 

lim f{R) oc R^. (79) 

R-^oo 

We see that the holographic /(i?)-gravity model (78) can satisfy the requirement of having to 
inflation (79). 

Replacing Eq. (76) into ECHDE density (36) one can obtain 



'^=4^"^' Vl2 



^ + Qfln 



(80) 



Solving the differential equation (75) for the energy density (80) yields the entropy-corrected 
holographic /(i?)-gravity model as 



2 



fiR) = XR^-c'R+[^yhln{R)-^ 



In 



\k^Rj 



2 



(81) 



which in the absence of correction terms, i.e. a = /3 = 0, recovers the result of holographic 
/(i?)-gravity model (78). As we already mentioned the correction terms in ECHDE density (80) 
become important in the early inflation era. Equation (81) also confirms that besides the term 
XR^, the corrections make sense during the inflation when R oo. 
The conformal time ij for the scale factor (73) yields 

= (82) 

Jo a aoH v / 

Here to obtain rj = ri(R) one cannot replace t by i? in (82). Therefore for the scale factor (73) 
one cannot obtain the /(i?)-gravity models corresponding to the NADE (47) and ECNADE (62) 
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densities. To avoid of this problem we set t ^ oo for the upper Umit of the integral (82). Hence 
the result yields 

-.2r>\ ^V2 



V 



f 

Jo 



dt 



(aoH) 



-1 



12 



Replacing Eq. (83) into NADE density (47) yields 



Substituting Eq. (84) in the differential equation (75) gives the solution 

f{R) = Ai?2 _ n'^alR, 



(83) 



(84) 



(85) 



where A is an integration constant. Note that the new agegraphic /(-R)-gravity model (85) like 
(78) satisfies the inflation condition (79). 

Replacing Eq. (83) into ECNADE density (62) one can get 



PA 



IP" 



R + 



ciqR 
~T2 



P + aln 



12 



k^alR 



(86) 



Solving the differential equation (75) for the energy density (86) yields the entropy-corrected 
new agegraphic /(i?)-gravity model as 



fiR) = XR'-n'alR+(^^ 



a 



f3ln{R)-- 



In 



12 



k'^alR 



(87) 



which for a = ^ = 0, the above result reduces to the new agegraphic /(i?)-gravity model (85). 
Here also like (81), during the inflation era not only the term XR^ but also the correction terms 
become considerable in Eq. (87). 



8 Conclusions 

Here, we considered the ordinary and entropy-corrected versions of the HDE and NADE models 
which are originated from some signiflcant features of quantum gravity. The HDE is motivated 
from the holographic hypothesis [4, 5, 6] and the NADE is originated form uncertainty relation 
of quantum mechanics together with the gravitational effect in GR [17, 18]. Among various 
candidates to explain cosmic accelerated expansion, only HDE and NADE models are based 
on the entropy-area relation. However, this definition can be modified from the inclusion of 
quantum effects, motivated from the LQG. Hence the ECHDE and ECNADE were introduced 
by addition of correction terms to the energy densities of HDE and NADE, respectively [15]. 

We investigated the HDE, ECHDE, NADE and ECNADE in the framework of /(i?)-gravity. 
Among other approaches related with a variety of DE models, a very promising approach to DE 
is related with the modified theories of gravity known as /(-R)-gravity, in which DE emerges 
from the modification of geometry. Modified gravity gives a natural unification of the early- 
time inflation and late-time acceleration thanks to different role of gravitational terms relevant 
at small and at large curvature and may naturally describe the transition from deceleration 
to acceleration in the cosmological dynamics [42]. We reconstructed the different theories of 
modified gravity based on the f{R) action in the spatially flat FRW universe for the three classes 
of scale factors containing i) a = ao{ts — t)~^, ii) a = aot^ and iii) a = aoe^* and according to 
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the original and entropy-corrected versions of the HDE and NADE scenarios. Furthermore, we 
obtained the EoS parameters of the corresponding /(i?)-gravity models. Our calculations show 

that for the first class of scale factor, the EoS parameter of the holographic and new agegraphic 
/(i?)-gravity models always crosses the phantom-divide line. Whereas for the second class, the 
EoS parameter of the mentioned models behaves like the quintessence EoS parameter. The EoS 
parameter of the entropy-corrected holographic and new agegraphic /(i?)-gravity models for the 
both of first and second classes of scale factors can accommodate the transition from quintessence 
state, ujR > —1, to the phantom regime, ujr < —1, as indicated by recent observations. For 
the third scale factor, the EoS parameter behaves like the cosmological constant. Also the 
/(i?)-gravity models corresponding to the HDE, ECHDE, NADE and ECNADE can predict the 
early-time inflation of the universe. 

Note that although f{R) theories offer a chance to explain the acceleration of the universe, 
they are not free of problems. As an alternative to DE for driving the late-time cosmic accel- 
eration, the /(i?)-gravity needs to pass the cosmological tests, including the constraints about 
the cosmic expansion and the cosmic structure formation. In addition, as a modified grav- 
ity theory, it needs to pass the solar system test, such as the constraints on the Brans-Dicke 
theory [43]. Besides, in many models of /(i?)-gravity, including the ones called realistic, the 
"fine-tuning" and the "cosmic coincidence" problems which are the two well-known difficulties 
of the cosmological constant problems, have not been essentially solved. Although there remain 
possibilities to solve the coincidence problem (see e.g. [44]), the existence of small parameters 
in /(i?)-gravity models is still the most important problem. In case of ACDM model, the scale 
of the cosmological constant is very small compared with the Planck scale, which is unnatural, 
especially from the viewpoint of the unified theory of the particle physics. Even in the realistic 
f{R) models like Hu-Sawicki's one [45] which satisfies both cosmological and solar system tests, 
there appear the small parameters, which would be still unnatural. 
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